On Free-Electron Laser Growing Modes and their Bandwidth 
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Free-electron lasers play an increasing role in science, from generating unique femtosecond X- 
ray pulses for single short recording of the protein structures to amplifying feeble interactions in 
advanced cooling systems for high-energy hadron colliders. While modern Free-electron laser codes 
can describe their amplification mechanism, a deep analytical understanding of the mechanism is 
of extreme importance for a number of applications. Mode competition, their growth rates and 
amplification bandwidth are among the most important parameters of a free-electron laser. A 
dispersion relation, which defines these important characteristics, can be solved analytically only 
for a very few simple cases. In this letter we show that for a typical bell-shape energy distribution 
in electron beam there is no more that one growing mode. We also derive an analytical expression 
which determines the bandwidth of the free-electron laser. 



Since the invention of free-electron lasers (FELs) by 
Madey the FEL field has made tremendous 

progress both experimentally ([|[, 3, @) and theoreti- 
cally. In an FEL electrons propagate through a wiggler, a 
device with a periodically oscillating transverse magnetic 
field, which causes the electrons to radiate at a resonant 
wavelength of 
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where a w = eA w /mc is the wiggler normalized vector 
potential, X w is the undulator period and 7 = £/mc 2 
is the electron relativistic factor. A resonant interaction 
between the electron beam and the TEM field of the ra- 
diation can become unstable, resulting in the exponential 
growth of the radiation power. This results in the gener- 
ation of coherent radiation. Such high-gain FELs are of 
great interests for current and future applications, and 
are the subject of the discussions in this letter. 

Since exact analytical solutions in three dimensions are 
typically intractable, numerical code such as GENESIS 
@ and GINGER Q are used for practical FEL designs. 
However, there do exist some analytical treatments such 

Q, and 

that can provide insight into the del- 
icate phase-space dynamics, statistical characteristics of 
FEL radiation, and can be used for testing the validity of 
the assumptions made in the FEL codes. The FEL linear 
regime describing the initial amplification of the density 
perturbation is frequently used for an in-depth analysis 
of the FEL physics, and most treatments assume a wide 
electron beam co-propagating with the TEM wave. 
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Still a number of fundamental questions about FELs 
remain unanswered. It is understood that the self- 
consistent ID FEL equation may have a variable number 
of modes depending upon the details of the energy distri- 
bution, with asymptotic formulae for the largest growing 
modes discussed in 11 [ and @ and a direct discussion of 
the diverging number of modes given in fl2l ] . For exam- 
ple, it is known that for a beam with a gaussian energy 
distribution an infinite number of modes exist. But there 
is no answer to how many of these modes are growing, 
or what frequency cutoffs might exist for these growing 
modes. 

Saldin, Schncidmiller and Yurkov provided the most 
general treatment of the high gain free-electron laser op- 
erating in the small-signal, linear regime [ij]. Using a 
Laplace transformation they reduced the self-consistent 
Maxwell- Vlasov equations [13[ to a dispersion relation of 
the form 
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where s is the Laplace transformation variable along the 
longitudinal direction, A 2 is a longitudinal space charge 

parameter, and D(s) is the dispersion integral given by 
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where we utilize the conventions in [14] of P = (£ — 
£o)/(p£o) as the normalized energy deviation, C = 
(k w — cj(1 + o4)/(2c7 2 ))-Lg the normalized detuning from 
the FEL resonance wave number, F the normalized en- 
ergy distribution of the electron bunch, the normalized 
to the e-fold gain length, Lq, and Pierce parameter, 
p = {k w L>G)~ l , and where k w = 2ir/\ w . The authors 
also found analytical solutions for the case of the mono- 
energetic and Lorentzian energy distributions as roots of 
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a cubic equation and expressed the solution for the sys- 
tem in the form 



(4) 



where the V n contain the initial conditions and the s n are 
the solutions of the dispersion relation. Most generally, 
for some energy distribution the solution is a linear su- 
perposition of exponential modes of the form e s,z where 
{St} are the solutions to the corresponding dispersion re- 
lation. 

Hence, the solutions with 5i(si) > is growing expo- 
nentially and knowing the number of these modes and 
their growth rate is the key for any high-gain FEL with 
gain much larger than unity. In this letter we will prove 
that for an electron beam with a typical bell-shaped en- 
ergy distribution, there is no more than one growing 
mode, and therefore in a high-gain FEL one can use the 
approximation that V(z) = Vbe s ° 2 . From the analyti- 
cal results for monocnergetic and Lorentzian beams, it is 
already known that there is only one growing mode. 

For most short wavelength FELs the longitudinal space 
charge does not play an important role and we neglect it 
in this letter, we therefore focus on the reduced dispersion 
relation 



w(s) = s- D(s) = 
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To study this, we use the Argument Principle of com- 
plex analysis. This is an analogy of the Nyquist stability 
criterion ([IH, [HI); the change in the argument of a 
complex function w{z) meromorphic on a closed contour 
C, is given by 



W = Z - P 
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where Z is the number of zeros of w{z) within the contour 
and P is the number of poles, counting multiplicity. 

Consider the example of the contour given in figure 
(|Tj), which encompasses the entire right half-plane in the 
R — > oo limit. By calculating the winding number, W, of 
the dispersion relation along this contour, we can make a 
general statement about the number of growing modes. 

First, we must determine the number of poles in the 
right half-plane of s. To determine this, look at the eval- 
uation of the dispersion integral D(s), since the poles 
of the full dispersion relation depend upon the poles of 
D(s). When evaluating the Laplace transform of the lin- 
ear equation for the FEL perturbation, it is assumed that 
5ft (s) > by consideration of causality. Specifically, 
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requires that, for t > 0, Jft(s) > for the integral to 
converge. When applied to the evaluation of D(s) by 
integrating over the complex P-plane, this means that 




FIG. 1: The contour over the complex-valued s that encom- 
passes the entire right half-plane. This contour determines 
the number of growing roots. 
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FIG. 2a: Contour integration for D(s) with K(s) > 0. 
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FIG. 2b: Contour integration for D(s) for the proper analytic 
continuation of 5R(s) over the entire complex s-plane. 
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the contour is given by figure (|2a|) . The correct analytic 
continuation is in figure (|2bj) as first shown in (l7j . 

The pole at is — C must remain above the contour 
integration when the dispersion relation is analytically 
continued to deal with the case of 3?(s) < 0. This causal- 
ity requirement means that the contour remains in the 
lower half-plane of P, and therefore the poles of D(s) are 
located at points given by the form 



<C + P* n ) 



(8) 



for a P* with an imaginary part less than zero. This 
translates to the poles of D(s) all existing in the left half- 
plane, so there are no poles inside our contour and P = 
0. Therefore the number of growing modes is given by 
the winding number of the dispersion relation! Causality 
in the Laplace transform is therefore directly related to 
the pole structure of the dispersion relation, which in 
turn gives us useful information about the number of 
amplifying modes. 

From very general considerations of the shape of F, we 
can arrive at a very solid schematic understanding of how 
the contour maps. The arc at infinity we parameterize 
by s = Re tB for 9 e (— n/2, tt/2]. Along this portion of 
the contour, D ~ Rr 2 and the dispersion relation is the 
identity map. 

Along the vertical end of the contour, s = it for t 6 
(oo, — oo ) and the dispersion integral is given by 



dP 



dF 



1 



V \ dP- 



.dF 



1 



dP i(t + C + P) 
- wF'(P = -t-C) 



(9) 



dP i{t + C + P) 
where we have used the well-known identity 

[■■■]-dx = V\ [•••]-] da + %Tt\->-\8(x)dx (10) 
x \ x J 

where V indicates the Cauchy Principal Value. The first 
term is purely imaginary since F and (t + C + P) ~ 1 are 
real functions, and the second term is pure real. The 
real part of the above expression is very simple, and pro- 
vides information when the mapping crosses the imag- 
inary axis, i.e. this happens when F' = 0. As a first 
example, we consider a simple bell curve shape for F, a 
curve with a single extremum F' = at t + C = t . 

In this case the map of the vertical line crosses the 
imaginary axis, $l(it — D(it)) = once and only once. 
Furthermore, because F has a bell shape, at the local 
maximum the second derivative of F - the derivative of 
the real part of the contour - is negative. Therefore, the 
contour is crossing from the right half-plane to the left 
half-plane. This leaves two possible options for the con- 
tour, given in figure (|3aj|3bp . The contour which crosses 
above the real axis corresponds to a winding number of 
one, while the contour which crosses below the real axis 
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FIG. 3a: Contour for a growing mode for a bell curve distri- 
bution. 
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FIG. 3b: Contour for zero growing modes. 



corresponds to a winding number of zero. Hence the con- 
tour in figure [3a] corresponds to a single growing mode, 
and the contour in figure[3blcorresponds to the case when 
all modes are either oscillating or decaying. Thus, we 
prove that for a typical bell-shaped energy distribution 
the one-dimensional FEL has one or no growing modes. 

Because the contour crosses the imaginary axis at 
t + C = to. the location along the imaginary axis, ei- 
ther above or below zero, depends upon the detuning. 
Thus, by inserting this solution into the imaginary part 
of the dispersion function, equation ([5]), we can write an 
expression for the detuning at which the growth rate of 
the FEL vanishes: 
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For a gaussian distribution given by 



F(P) = 




-P 2 /(2a P )\ (12) 



this gives the critical detuning as 

& = -\ (13) 

Up 

More generally, it is clear by dimensional analysis that 
for any energy distribution with only one scale ap the 
critical frequency scales as C* oc o-p 2 with the numeri- 
cal quantity on the order of unity and dependent on the 
particulars of F'(P). 

To summarize this result, assuming only that F(P) 
has a single maximum and goes to zero at P — > ±oo, we 
have shown that causality demands a single mode which 
amplifies frequencies with C > C* oc —o-p 2 with some 
numerical quantity depending on the specifics of the dis- 
tribution and ap is the energy spread parameter for the 
distribution. We should note that this result holds even 
for the cold beam with zero energy spread, which am- 
plifies all frequencies above the FEL resonance frequency 
(note that negative C corresponds to higher frequencies). 

It is worth mentioning that, generally speaking, an 
FEL with an energy distribution with N humps can have 
up to N growing modes. It also means that for any typi- 
cal energy distribution there is a finite number of growing 
modes. The arguments used to construct the contours 
in figures (|3aH3b| may also be used to draw schemat- 
ics for a multiple-peak energy distribution. Every local 
maximum represents the contour crossing from the right 
half-plane to the left half-plane and every local minimum 
represents the contour crossing from the left to the right 
half-plane. Extrema that are neither maxima nor min- 
ima (i.e. plateaus) are bounces and do not contribute. By 
looking at the sign flips in the imaginary part, analogous 
diagrams can be drawn and by a similar argument there 
are zero, one, or two growing modes corresponding to one 



mode for each resonant frequency that corresponds to a 
local maximum, and having one, both, or none of them 
being amplified at a given frequency. Similar criterion to 
the above for a critical detuning may be developed us- 
ing the same arguments. In practice this is of little real 
interest, but it is intellectually satisfying that this formal- 
ism can consider multiple resonant frequencies caused by 
multiple peaks in the energy distribution. Because the 
specifics of the number of growing modes is independent 
of the energy distribution beyond some basic considera- 
tions on the shape of the distribution, we argue that this 
result may be regarded as topological. 

The method used to obtain this result does not depend 
on the fact that we are dealing with an FEL, and indeed 
these results have been obtained in frequency space by 
Nyquist (lo| . Penrose for a stationary plasma driven by 
an external electric field with frequency oj (l8j . and for 
general "growing waves" or other plasma instabilities 
[20j . This particular method appears to be periodically 
rediscovered by physicists over the decades. Because the 
Laplace transform and Fourier transform are related by 
a Wick rotation, the same formal structure exists either 
looking at the initial value problem or at frequency re- 
sponses. 

Therefore we have proven that, for a reasonable bell- 
shaped type energy distribution in the electron beam, 
the FEL dispersion relation has only one growing solu- 
tion, which exists only within a frequency range defined 
by a simple expression (jlip . We extend our qualitative 
topological considerations for energy distributions having 
multiple peaks and conclude that for a smooth and finite 
energy distribution there can only be a finite number of 
growing modes, simply related to the number of extrema 
of the energy distribution. 
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